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Can the dark matter annihilation signal be significantly boosted by substructures? 
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A very general cosmological consideration suggests that, along with galactic dark matter halos, 
much smaller dark matter structures may exist. These structures are usually called ’clumps’, and 
their mass extends to 10~^Mq or even lower. The clumps should give the main contribution into the 
signal of dark matter annihilation, provided that they have survived until the present time. Recent 
observations favor a cored profile for low-mass astrophysical halos. We consider cored clumps and 
show that they are significantly less firm than the standard NFW ones. In contrast to the standard 
scenario, the cored clumps should have been completely destroyed inside ~ 20 kpc from the Milky 
Way center. The dwarf spheroidals should not contain any dark matter clumps. On the other hand, 
even under the most pessimistic assumption about the clump structure, the clumps should have 
survived in the Milky Way at a distance exceeding 50 kpc from the center, as well as in low-density 
cosmic structures. There they significantly boost the dark matter annihilation. We show that at 
least 70% of the clumps endured the primordial structure formation should still exist untouched in 
the present-day Universe. 

PACS numbers: 95.35.-|-d; 95.30.Cq; 98.52.Wz 


I. INTRODUCTION 

A very general cosmological consideration suggests 
that, along with large dark matter haloes, corresponding 
to visible astronomical objects, like galaxies or galaxy 
clusters, much smaller dark matter structures may ex¬ 
ist. Indeed, the inflation inevitably leads to a Zel’dovich- 
Harrison-like spectrum of the primordial cosmological 
perturbations, which is flat and contains perturbations 
of any scale with approximately the same amplitude [I]. 
The perturbations with masses lower than some cutoff 
value Mcut may be destroyed in the early Universe by 
the free streaming. If the dark matter (hereafter DM) 
is cold and non-interacting, which is the most popu¬ 
lar hypothesis now, Mcut is mainly determined by the 
mass of dark matter particle rridm- Roughly speaking, 
Mcut where Mpi is the Planck mass. For 

'nidm = 100 GeV, different estimations of Mcut vary from 
1O“^^M0 [2] to 1O“®M0 [3]. We denote the clump virial 
mass and radius by md , rd ■ 

We consider the DM particles (hereafter DMPs) of 
mass 100 GeV not by accident. If we suppose that dark 
matter was in thermodynamical equilibrium in the early 
Universe, current DM abundance allows estimating of the 
DM annihilation cross section {<Jv)a — 3 • 10“^®cm^/s. 
This is a typical value for weak interactions, which im¬ 
plies that DMPs may interact weakly. Of course, other 
explanations of the dark matter abundance are possible, 
and the ’correct’ cross section value can be just a coinci¬ 
dence. However, the weakly-interacting massive particles 
(WIMPs) are one of the most popular dark matter candi¬ 
date now. The WIMP mass in realistic models typically 


exceeds 100 GeV. A lighter WIMP is possible, but a theo¬ 
rist should take special precautions to avoid tensions with 
the LHG results, since a light WIMP is usually accompa¬ 
nied by respectively light charged particles in theories like 
supersymmetry. Thus the existence of tiny dark matter 
halos (traditionally named ’clumps’) is not a theoretical 
assumption, but an inevitable property of many models 
of the cold dark matter, if DMP mass exceeds ^ 1 MeV. 

The correct estimation of the dark matter dumpiness 
is of critical importance, especially for the indirect dark 
matter detection. Indeed, the number of annihilations in 
a volume dV in an interval of time dt is [32] 

hav)^dTdV ( 1 ) 

2 ml 

where (cv) is the averaged multiplication of the DM an¬ 
nihilation cross-section on the relative velocity and is 
the DMP mass. Thus the annihilation signal is propor¬ 
tional to the so-called boost factor C = (p^)/(p)^, which 
takes into account the enhancement of the annihilation 
signal by dark-matter substructures. In the standard cos¬ 
mological scenarios the smallest halos collapse the first, 
and therefore have the highest density contrast, mak¬ 
ing the main contribution to the boost. We can define 
a characteristic radius of each halo where its density 
profile gets steeper than r“^, i.e., dlogp/dlogr = —2 
at Ts- Then we may introduce the halo concentration 
Cuir = Ts/Rvir- In the case of NFW halo, this definition 
coincides with the standard one. 

The information about the smallest DM structures is 
poor at recent. First of all, astronomical observations 
of the clumps are impossible, since their gravitational 
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potential is too small to hold baryon matter. The the¬ 
oretical view on the clump structure and velocity distri¬ 
bution of the DM particles is also rather uncertain. The 
clump mass is very small (< 1 O“®M 0 for neutralino dark 
matter) comparing to the galaxy masses, and results and 
relationships obtained in standard N-body simulations 
should be extrapolated on many orders of magnitude [4]. 
Moreover, as we will see below, the present-day number 
of clumps strongly depends on their interaction with var¬ 
ious baryon objects, like stars, galaxy disc etc. It makes 
a reliable simulation of clump abundance very difficult. 

The main fraction of the clumps formed in the early 
Universe has not survived to our time. There are many 
mechanisms of clump demolition, and they can all be di¬ 
vided into two groups. A significant part of clumps was 
tidally destroyed almost immediately after their occur¬ 
rence by larger halos during the hierarchical structure 
formation. We will imply the standard ACDM model 
hereafter in this paper. In the framework of this scenario, 
the hierarchical clump destruction is not sensitive to the 
clump density profile. Indeed, the spectrum of primor¬ 
dial perturbations is flat in this case, and objects of dif¬ 
ferent masses collapse almost simultaneously. The time 
interval between the occurrence of the smallest clumps 
and more massive objects that capture the clumps and 
destroy them tidally is smaller than the clump viriliza¬ 
tion time. Consequently, even if present-day clumps have 
cusps, they were formed later. Therefore, we may use a 
simplified estimation of the hierarchical clump destruc¬ 
tion by calculating an energy gain per each tidal inter¬ 
action and the number of tidal interactions [5]: a clump 
is considered to be destroyed if its internal energy in¬ 
crease due to the tidal shocks exceeds its total energy 
\E\ ~ GM/R. This estimation leads to the oc dM/M"^ 
mass spectrum of the clumps survived the structure for¬ 
mation and suggests that only about 0.1 —0.5% of clumps 
endure the hierarchical tides in the each logarithmic mass 
interval. N-body simulations [ 6 ] support the mass spec¬ 
trum, but suggest a several times higher clump survival 
ratio. However, even 0.1 — 0.5% of clumps endured are 
sufficient to provide a very significant boosting [7]. 

The second large group of clump destruction mecha¬ 
nisms can be named astrophysical. The clumps can be 
disrupted by separate stars, galaxy stellar discs, tidal per¬ 
turbations in the galactic gravitational field, gas clouds 
etc. Contrary to the hierarchical destruction, this pro¬ 
cess goes on up to now, and its efficiency is critically 
defined by the internal structure of the clumps. We will 
consider these processes in the present work and try to 
show that the clumps may be much less persistent than 
it is generally believed. 


II. ENERGY STRUCTURE OF CLUMPS 

It might appear at first sight that the difference be¬ 
tween a cuspy and a cored profile is not that large, es¬ 
pecially if the core is small. However, this conclusion is 


incorrect, if we take into account a possible velocity dis¬ 
tribution anisotropy of DMPs. In order to illustrate it, we 
consider the evolution of DMP energy distribution dur¬ 
ing the halo formation. The halo is formed from a linear 

initial perturbation. The potential well of it was shallow. 
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As a result, the initial total energies e = the par¬ 

ticles forming the halo were close: under very general as¬ 
sumptions, they lay in the interval e € [—1.5$^^^.; —^vir] 


(where 

itvir 


, i.e., they differed no more than 


1.5 times [ 8 ]. The final potential well is much deeper. 
The ratio of its depth (i.e., the central potential (j){0)) to 
the potential on the virial radius —^vir is (j){0)/^yir ~ 30 
for the Milky Way [9] and even larger for smaller objects, 
since they should be more concentrated. 

It is convenient to characterize each particle, instead 
of its energy, by its apocenter distance ro, which is the 
largest distance that the particle can move away from 
the center, ro is defined by the particle energy and the 
shape of the potential well. For instance, if the par¬ 
ticles had kept their initial energies from the interval 
[-1.5<I>„r; -‘bmr], their ro would lie in tq G 

A question appears: how strong can the energy evolu¬ 
tion of the system during its formation be, i.e., can the 
relaxation be arbitrarily strong, or the ratio between the 
final e/ and the initial Cj energies of the particles is some¬ 
how limited for the majority of DMPs? In practice, two 
scenarios are possible: 


1. Strong energy relaxation. The particle energies 
are completely redistributed, final particle energies 
cower all the potential well [(()(0); —Then 
ro of the particles cover all the possible interval 
[ 0 , Rvir] more or less evenly. 

2. Moderate energy relaxation [ 8 ]. The final total spe¬ 
cific energy e/ of most of the particles differs from 
the initial ones no more than by a factor k 

— < A: (2) 

e* 

and k <C (j){0)/^vir- There can be particles that 
have changed their energy stronger, but their frac¬ 
tion is respectively small. 

The relaxation is moderate if fc < Cvir/4: [10]. This 
condition is very strict for galaxy clusters, but quite soft 
for low-mass galaxies with Cyir > 20 (though the real 
density profile may significantly differ from the NFW one, 
we used here the NFW halo concentration Cyir because of 
its popularity and in view of the fact that characteristic 
values of Cmr for various types of astronomical objects 
are well known). 

The moderate relaxation leads to a very important 
consequence. Since initial particle energies were close 
to —^vir and changed no more than k times during the 
halo formation, the halo contains only a small fraction 
of particles with the energies lower than —k^yir (it is 
important that \k^yir\ <C |(()( 0 )|), which roughly corre¬ 
sponds to ro — Ryir/k. Particles with more compact 
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orbits are almost absent in the halo. It means that the 
density profile in the center of the halo, inside Ryir/k, 
is mainly formed by the particles that arrive there from 
the outside. It leads to a very peculiar structure of the 
clump, differing drastically from the clumpy one. 

First of all, it is easy to show that a cusp may form 
only in scenario 1, i.e., if the relaxation is strong [ 8 , 10]. If 
the relaxation is moderate, the density profile obligatory 
has a core. Nevertheless, the profile still can be very 
concentrated: the ratio of the core radius to the virial 
one can be arbitrarily small. The velocity distribution in 
the case of moderate relaxation is isotropic only in the 
very center of the halo [ 11 ] and gains some anisotropy 
with radius. 

What scenario, 1 or 2, is realized in nature? As of now, 
the answer to this question is not known with certainty. 
The mass hierarchy of the real Universe is so vast (from 
rich galaxy clusters of mass ~ IO^^Mq down to clumps 
of mass 1 O“®M 0 or even lower) that it can be covered 
in no simulation, and a far interpolation is necessary. It 
makes discussible some simulation results concerning the 
clumps. For instance, if we interpolate the power-law 
grows of halo concentration Cvir with its mass decreasing 

[ 12 ] , we obtain extremely large concentrations for clumps. 
This result is discussable [4]. Recent N-body simulations 

[13] suggest that the power-law grows of Cyir breaks for 
small-mass halos. Moreover, the simulations may suf¬ 
fer from underestimated numerical effects in halo centers 

[14] . N-body simulations favor strong energy relaxation 
leading to the cuspy profile. Their results nearly approx¬ 
imate galaxy clusters. However, if we consider smaller 
halos, like spiral galaxies or dwarf spheroidals, the model 
of moderate relaxation describes observations much bet¬ 
ter [ 8 , 10]. In particular, observations suggest cored pro¬ 
files [15-20]. Therefore, we may expect that formation of 
low-massive clumps also follows the moderate relaxation 
scenario. 

The moderate energy evolution can be characterized by 
the above-mentioned relaxation parameter k. As we will 
see in the Comparison of the predictions of core and cusp 
models section, the larger k is, the stronger the relaxation 
is, and the more stable the clumps are. The clumps are 
the most fragile if fc = 1 , i.e., when the clumps are formed 
with a negligible relaxation, and therefore their particles 
have To ~ rd- It is reasonable to expect that the moder¬ 
ate relaxation is much stronger (k cs. Cyirjdt) [10]. How¬ 
ever, in this paper we will mainly consider the case k = 1 
for simplicity, as it gives the most pessimistic estimation 
of clump survival. Thus, though the choice between sce¬ 
narios 1 and 2 is now unclear, the model we consider in 
this paper may be used to estimate the minimal value of 
the boost factor. 

Moreover, if the assumption of the moderate energy 
evolution is true, it suggests that the majority of the 
particles still concentrate near e = the energy 

exchange of the particles during the relaxation is more 
or less a stochastic process, while the total energy of the 
system should conserve. Thus the narrowness of the ini- 



FIG. 1: The profile of a clump formed in the case if its energy 
relaxation was moderate (solid line). The NFW profile with 
the same concentration parameter Cvir = 25 is reproduced 
for comparison (dashed line). The density is represented in 
arbitrary units. 


tial energy distribution, the relative smallness of the en¬ 
ergy evolution, and the fact that the final potential well 
is much deeper, than the initial one, together result in 
crowding of the particle apocentre distances rg near Rmr , 
i.e., (ro) ~ rd. 


III. TYPICAL CLUMP PARAMETERS 


If the dark matter particle is a WIMP of mass ex¬ 
ceeding 100 GeV, the minimal clump mass is roughly 
equal to 1 O“®M 0 , which corresponds to Ryir ~ 0.05 pc. 
The average density p of the clump is defined from 
p ■ = Myir. We will use the following des¬ 

ignations: mg = Myir/lt)~^MQ, rg = Ryir/O.Q^ pc, 
Pq = 1.9 • 10“^ Mq/pc^. For a clump with mass 
Myir = 1 O“®M 0 we have mg = rg = pg = 1 . 

As we could see, almost all the particles in our model 
have ro Ryir, and the dynamical time of the clump 

can be estimated as tf ~ l/If 

y Cr^vir 

Myir = 10 ~^Mq, tf ^ 10® years. In the standard 
ACDM cosmology p depends on the clump mass only 
logarithmically [ 1 , 10 ], and as a result, the clumps of all 
masses have comparable dynamical times. The time is 
significantly larger than typical times of any clump de¬ 
struction mechanisms. Therefore, the clump disintegra¬ 
tion can be considered in the instant tide approximation. 
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IV. CLUMP DESTRUCTION: THE CASE OF 
THE MOST FRAGILE CLUMPS k = 1 

There are several mechanisms of clump destruction in 
our Galaxy. We will consider only the most important 
ones. As we showed at the close of section Energy struc¬ 
ture of clumps, the case k = 1 corresponds to ro ~ rd and 
gives the most pessimistic estimation of clump survival. 
We will actively use this fact in the current section. 


A. Destruction on halo stars 


The clumps can be destroyed by collisions with stars. 
Let us denote the average star mass, star number density 
and star velocity dispersion by m*, n*, and cr*, respec¬ 
tively. The dispersio n of clum p velocities is a^- We 
may introduce cr = a/ct* -I- . The characteristic time 

of clump destruction by collisions with stars is [21, eqn. 
8.55] 


^ _ 17 amcirh 
200 Gmln^^rQ ’ 

Since ro ~ rd , we obtain: 

^ 17 amd 
200 Gmln^^rd ’ 


(3) 

(4) 


The main fraction of time the clumps spend in the Milky 
Way halo. The star concentration there can be estimated 
[7] as n* = (p?i/m*)(r©/r)^, where r© = 8.5 kpc, ph = 
1.4 • lO'^Mg/pc^, m* ~ 0.4M©. We should integrate n* 
over the clump trajectory, from the minimum r^m to the 
maximum r^nax radius of its orbit. However, n* rapidly 
drops with radius, and the region near r^in makes the 
main contribution into the clump destruction. Assuming 
a = 200 km/s, we obtain from (4): 

td ^ ■ 8 • 10^°years. (5) 


As we can see, td exceeds the age of the Universe even 
for low-massive objects. Only very small clumps Myir < 
10“”^M© can be destroyed by the halo star collisions. 


with Md = 8 • lO^'^M©, Td = 4.5 kpc. Dividing the total 


clump energy \E\ 


Gwid 


rd 


by (6), we obtain the number 


of passing Nd needed to destroy the clump by the disc 

( 8 ) 


_ 3r>i ^ 
SGMl rl 


Substituting[33] here n_L ^ 100 km/s, we obtain Nd — 
4.5 • . As we can see, Nd depends only on 

the average density of the clump. Since p depends on the 
clump mass only logarithmically in the standard ACDM 
model [1, 10], the clump destruction by the disc depends 
almost not at all on the clump mass. 

The motion of a clump in the Galaxy can be charac¬ 
terized by the minimum rmin and the maximum rmax 
radii of its orbit. Equation (8) shows that a single pass¬ 
ing through the disk is enough to disintegrate a clump 
if fmin < 17 kpc. So high efficiency of the mechanism 
under consideration in the central area of the Milky Way 
justihes our neglect of other clump destruction factors, 
like the galactic bulge, gas nebulas etc. In the regions 
where they are effective, the clumps are disrupted any¬ 
way. 

On the other hand, the disc surface density rapidly 
drops with radius, and so does the efficiency of clump de¬ 
struction. A clump passes through the disc twice in each 
revolution. Let us consider a clump with a circular orbit 
of radius 25 kpc. The orbital period of the Solar System 
about Galactic Genter is Tmw ~ 230 Myr. If we accept 
that the Milky Way density profile is p oc (which 
is a very rough approximation), T oc r [22], i.e., the or¬ 
bital period of the clump is T ~ 230Myr ■ 25kpc/8kpc ~ 
700 Myr. The clump passed through the disc ~ 30 
times in the lifetime of the Galaxy. According to Eqn. 8, 
the clump should be destroyed in ~ 30 passages. Thus 
r = 25 kpc is approximately the border of the region 
where the clumps were disrupted by the disc in the life¬ 
time of the Galaxy. Clumps with larger r also can be de¬ 
stroyed, if they have an elongate orbit and pass through 
the disc close to the center. However, the fraction of these 
clumps falls approximately as and becomes negligible 
at r ~ 50 kpc. 


B. Destruction on the Milky Way disc 


However, the clumps can be disrupted by the tidal per¬ 
turbations of the galactic disc as a whole. A single pass¬ 
ing of a clump through the disc changes its energy on 


6E = 


IbTT^G^o^rg 


3v 


2 

± 


( 6 ) 


where Od is the surface density of the Milky Way disc, n_L 
is the perpendicular to the disc component of the clump 
velocity. We assume [7] 


Od{r) 


Md 

27rr2 




(7) 


C. Destruction by the gravitational field of the 
Galaxy 


If a clump has a non-circular orbit, it periodically ap¬ 
proaches the galaxy center and expresses tidal shocks. A 
single shock increases the clump energy on [21, eqn. 8.43] 


8E = 


4G^M^ir„ 

3u2ri,v. 




(9) 


where v is the clump velocity at rmin, and Mpijmin) 
is the total mass of the perturber (in our case, of the 
Milky Way) inside rmin, and (r^) is the mass-weighted 
mean-square radius of the clump. If the clump density 
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profile is p oc r~^ and rg = rd, (r^) = rg/3 = 
Though a real clumps apparently have other profiles, we 
may use approximation (r^) = r^;/3 in our estimative 
calculations. Dividing the clump energy E cii Gnici/rci 
by 5E^ we obtain the number of rotations necessary to 
destroy the clump 


that we consider the most pessimistic scenario). A gray 
zone lies between 20 and 50 kpc, some clumps should be 
destroyed there, some may survive. No mechanism ef¬ 
fectively disrupts clumps out of 50 kpc from the Milky 
Way center, even under the most pessimistic assumptions 
about the inner clump structure that we consider. 


Nd 


9 / "C Vfnin \ / '^cl \ 

4 [gmJ 

^ rain )) \ ^cl ) 


( 10 ) 


D. Dwarf spheroidals. Segue 1. 


If the orbit has a significant eccentricity, we may estimate 

"C — GIVlp{T<Yyt,iYi^ / 


NAr) = 2 ^. , 11 ) 

where Pp{r) is the average density if the perturber inside 
radius r. Thus the clump can be effectively destroyed if it 
passes through the areas where its average density is com¬ 
parable with the density of the perturber. The destruc¬ 
tion process again depends only on pd and therefore is 
almost insensitive to the clump mass. According to equa¬ 
tion ( 11 ), a clump should have orbit radius r ~ 18 kpc 
do be destroyed in one turn about the Galaxy center and 
r ~ 50 kpc to be destroyed in the lifetime of the Galaxy, 
for the conventional Milky Way models. 

We should underline that the mechanism efficiency 
strongly depends on the orbit eccentricity: a clump on a 
circular or near-circular orbit expresses no shocks. How¬ 
ever, there is a clump destruction mechanism that is not 
sensitive at all to the shape of the clump orbit. Indeed, 
let us consider a particles belonging to clump that ro¬ 
tates around its center having the apocenter distance rg. 
The tidal force from the Galaxy gravitational field can 
be estimated as {GMpir)/r^) ■ (rg/r) in the framework of 
the clump center, where r is the clump distance from the 
Galaxy center. The work of the force on the particle or¬ 
bit is ~ rg • {GMp{r)/r‘^) ■ (rg/r). If it is comparable with 
the bind energy of the particle ~ Gmd/ro, the particle 
can be teared out of the clump. We assume rg ~ rd, as 
we did before. The condition of clump destruction is 


Mpjr) 

p3 


rud 



or pd < Pp{r). 


( 12 ) 


If we compare equation (11) with the last one, we can 
see, that it is approximately equal to the condition for 
a clump to be destroyed in one or two turns around the 
Galaxy center. Thus the direct tidal disruption is in gen¬ 
eral less effective than the tidal shocks. On the other 
hand, equation ( 12 ) guarantees that all the clumps (in the 
most pessimistic case for their survival, when rg ^ rd) 
are destroyed inside r ^ 20 kpc, disregarding their orbits. 

To summarize: The clumps in the Milky Way halo 
are destroyed by the disc and by tidal effects in the 
Galaxy gravitational field. The efficiency of these two 
mechanisms is comparable. Collisions with halo stars 
are not effective. All the clumps should be destroyed in¬ 
side ~ 20 kpc from the Galaxy center (we should remind 


Among all astrophysical objects containing dark mat¬ 
ter, dwarf satellites of the Local Group draw special at¬ 
tention. They are very promising for detecting dark- 
matter annihilation. Though the dark-matter density 
in these objects is significantly lower than that at the 
Galactic center, they typically contain no sources of cos¬ 
mic rays, and therefore almost any high-energy signal 
detected from the dwarf galaxies could be the signal of 
dark-matter annihilation. 

Segue 1 draws the most attention among numerous 
satellites in the Local Group [23]. It combines several 
advantages: it is the closest satellite to Earth, it is lo¬ 
cated at high Galactic latitude (and, consequently, su¬ 
perposes on a low gamma-ray background) and expected 
to produce the strongest dark-matter signal [24], though 
the last statement is still discussable. Recent observa¬ 
tions can essentially constrain the most probable WIMP 
dark matter models (with the cross-section correspond¬ 
ing to the dark-matter abundance (ctu)o — 3-10“^®cm^/s 
and DMP mass range that has not been excluded by 
LHC experiments) only if the dwarf spheroidals are ex¬ 
pected to have a significant boost factor [25]. Therefore 
the question of the substructures in dwarf spheroidals is 
extremely important. 

Segue 1 has the highest known mass-to-light ratio 
of any observed galaxy, which implies that Segue 1 is 
strongly dominated by dark matter. Its density distribu¬ 
tion can be modelled by the Einasto profile 



with pe = 0.11 M 0 pc“^, Te = 150 pc, and n = 3.3 [26]. 
Unfortunately, all the parameters are poorly known and 
contain huge uncertainties: all the stars (~ 66 ) observed 
in Segue 1 lie inside 10' (~ 67 pc) from the center, which 
is two times smaller than rg. For instance, a completely 
different set of parameters pe = 0.094 Mq pc“^, rg = 
70 pc fits the observations equally well [27]. 

Apparently, Segue 1 does not contain a disc. However, 
the clumps could still be disrupted by stars or by the 
tides in the gravitational field of the dwarf. We consider 
both these mechanisms. 

Segue 1 contains approximately a thousand stars [28] 
that are situated inside 67 pc from the center. We may 
very roughly estimate n* ~ 8 • lO^'^pc"^. Substituting 
TO* = O. 4 M 0 and cr = 60 km/s (which is approximately 
the escape velocity from the center of Segue 1 [25]), we 
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obtain 

td — ■ lO^years. (14) 

As we can see, the smallest clumps are destroyed by the 
stars. However, clumps of mass me > 30, i.e., m > 
3 • 10“®M©, cannot be disrupted by this mechanism. 

In contrast, the clump destruction by the gravitational 
field of Segue 1 is extremely effective. We can see from 
eqn. (13) that the equality pd = Ppir) occurs at the 
distance ~ 1300 pc from the center of Segue 1 (if we 
suppose that mass distribution (13) is still valid there). 


V. THE CASE OF ROBUST CLUMPS k = c„i^/4 


The case of /c = 1 is extreme and not very probable, 
even if the relaxation is moderate. One may expect that 
k depends on Cyir of the halo. The most probable value 
of k is Cyir/4: [10], which is close to the maximum value of 
k, at which the relaxation may still be named moderate. 
Therefore we should consider the case k = Cyir/i. 

First of all, now k depends on Cyir- If fc = Cyir/4, 
(ro) > Tci/k = 4Ryir/cyir [8], and we should substitute 
this value to (9) and (6). Recent simulation suggest that 
the power-law growth of Cyir with the halo mass decreas¬ 
ing breaks, and the concentration of low-mass objects 
grows much slower. A sophisticated analysis [13] has 
shown that even the smallest halos have Cyir — 60. We 
will use this value as the most plausible. 

Equation (8) for the number of rotations necessary to 
destroy the clump by the Milky Way disc transforms into 


128 GMI 


(15) 


A reasoning similar to that used below equation (8) shows 
that r = 15 kpc is approximately the border of the re¬ 
gion where the clumps were disrupted by the disc in the 
lifetime of the Galaxy in the case of robust cored clumps 
with k = Cyir j4. This value is substantially smaller than 
that for the most fragile clumps fc = 1, but larger then the 
value obtained for the cuspy clumps r ~ 8 kpc [7]. Thus 
the durability of the cored clumps may vary in a wide 
range. However, an essential difference between cuspy 
and cored models remains. 

Equation (12) transforms into 


< Pp{T). (16) 

This equation shows that the dwarf satellites like Segue 
1 are so dense that it is difficult to expect that any cored 
clumps could survive in it. Indeed, if we accept Cyir — 60 
for clumps, then the clumps are destroyed in the region 
where /jp(r) > 225pd- Equation (13) shows that all 
the clumps should be destroyed inside ~ 2re, i.e., inside 
300 pc from the center of Segue 1. 

The conclusion of absence of any DM substructures is 
probably valid for all dwarf satellites of the Local Group. 
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FIG. 2: The annihilation signal J from the Milky Way as a 
function of the angle 9 between the observation line and the 
direction towards the Galaxy center. The normalization of J 
is arbitrary. The solid line corresponds to the case when the 
clumpy structure disappears at rdeatr = 50 kpc, the dashed 
line - to rdestr = 10 kpc. The boost factor at 50 kpc is set to 
be equal to C 50 = 20. 


The satellites, as we observe them now, are the most 
probably the central and the densest parts of ancient 
dwarf galaxies tidally destroyed by the Milky Way and 
M31. As a result, their density is very high, and the 
clumps are easily destroyed when approaching their cen¬ 
ters. We should also take into account that the tidal per¬ 
turbations destroying the satellites disrupted their clump 
structure as well. 


VI. COMPARISON OF THE PREDICTIONS OF 
CORE AND CUSP MODELS 

Now we can compare the survival of the cored clumps 
formed under the assumption of moderate relaxation 
(hereafter MR-clumps) with one of the standard model 
NEW clumps [5, 7]. As we can see, clumps are destroyed 
much easier in our case. 

This result should not surprise. The clumps of the 
standard model are significantly more stable not be¬ 
cause of the existence of cusp per se: MR-clumps can 
also be very concentrated and have large central density 
(Eig (1)). The principle difference is the particle orbits. 
The particles forming the cusp in the NEW model are 
confined in the cusp: they have small quasi-circular or¬ 
bits with the apocenter distance tq cs 0. According to 
eqn. 6, the energy perturbations are proportional to Tq. 
This fact along provides the cusp stability. Moreover, 
the particles with ro ~ 0 have a very short dynamical 
time, the tidal perturbations are adiabatic for these parti¬ 
cles, and equations like (8) or (11) are not valid anymore. 
The adiabatic invariant conservation makes the cusp ex¬ 
tremely resistant to tidal perturbations. A NEW clump 
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easily looses its outer layers, but the cusps are almost 
immune. They can be destroyed only in the bulge of the 
Milky Way, at the distances less then ~ 3 kpc from the 
Galaxy center [7], and most likely survive in the dwarf 
spheroidals. Since it is cusps that makes the main con¬ 
tribution into the boost factor, it should be significant in 
the standard model. 

The MR-clumps are built completely differently. They 
still can have a very high central density, and their pro¬ 
file can be relatively close to the NFW, except of a small 
central region, where they have cores (see Fig (1)). How¬ 
ever, the majority of the DMPs have the apocenter radius 
To comparable with the virial radius of the clump. It is 
true even for the core itself: The bulk of the particles 
only ’visit’ the halo center, having tq ^ Tc, where Tc is 
the core radius. As we could see, the particles with large 
apocenter distance tq are much easier for the tidal effects 
to detach. The MR-clumps have a much larger apocenter 
radius (ro) and therefore are much easier to be destroyed. 

Evidently the largest possible value (ro) = Rvir, which 
corresponds to the relaxation parameter k = 1. It is 
exactly the assumption we used in obtaining of the results 
of this paper. Therefore we have obtained the low limit 
on the clump stability. The most probable value of k is 
k cs Cmrl4:, where Cvir is the NFW halo concentration 
[10]. Then (ro) ~ 2Ryir/cmr, and the clump survival is 
somewhat higher. However, it is still much lower for the 
cored clumps than for the cuspy ones. 

We have shown that the clump survival strongly de¬ 
pends on their internal structure, which is not quite clear 
now. Since the annihilation signal rapidly increases with 
the DM density growth, compact and dense DM regions 
attract the main attention of the indirect DM searches. 
Our calculations show that the clump structure in these 
objects can be utterly destroyed under quite natural cos¬ 
mological assumptions. 

It is important that, even under the most pessimistic 
suppositions about their structure, the clumps in the 
Milky Way should survive at the distances exceeding 
50 kpc from the Galaxy center. It means that the clumpy 
structure should survive in almost all the dark matter in 
the Universe. If we consider the Local Group, ~ 80% 
of its DM containment form the halos of the Milky Way 
and M31 outside 50 kpc from their centers and proba¬ 
bly a common DM envelope of the system [29]. Broadly 
speaking, the main part of the dark matter in the Uni¬ 
verse forms low-density structures, like DM envelopes of 
galaxy clusters, walls and filaments of the large scale 
structures etc. All this dark matter should have kept 
entirely its clumpy structure. This fact is very impor¬ 
tant for the estimations of the extragalactic background, 
possibly created by the DM annihilation. 


VII. PHENOMENOLOGY: PREDICTIONS OF 
THE ANNIHILATION SIGNAL AND 
UNCERTAINTIES 

We should underline that the phenomenological conse¬ 
quence of the DM substructures is not only the annihi¬ 
lation signal boosting by some factor. The signal depen¬ 
dence on the DM density p changes. If the dark matter is 
homogeneous, the annihilation signal is proportional to 

. If the annihilation in clumps dominates, the signal 
is proportional to the number density of the clumps, i.e., 
to p. Thus, if the clumps exist in the present-day DM 
structures, the DM annihilation should be boosted, and 
the boost factor should be roughly proportional to p~^. 

Let us discuss the uncertainties and possible observ¬ 
able consequences of the model. As we have shown, the 
clumps survive outside of 50 kpc from the galactic center 
even under the most pessimistic supposition about their 
structure. On the other hand, [7] showed that even rather 
robust NFW clumps are strongly destroyed closer than 
8 kpc from the center. We may conclude that 20 — 50 kpc 
is quite a reliable estimation of the maximum distance 
from the Milky Way center where the cored clumps can 
still be destroyed. The boost factor estimations are much 
more ambiguous. 


A. A sketch of the standard theory of the boost 
factor 

If the comoving number density of halos in the mass 
range of [M; M -|- dM] at the redshift z is dp{M, z)/dM, 
the boost factor of the system C can be represented as a 
multiplication (see details in [30]): 

C = A-F,M{z)-[Bhalo]. (17) 

Here A is the mean halo overdensity defined by the cos¬ 
mological model, of a halo in units of the cosmic mean 
matter density, Fcoii(z) = f -^MdM/(p) is the mass 
fraction collapsed into cosmological halos at z {(p) is 
the average density of the Universe). Bhaio is a sin¬ 
gle halo clumping factor Bhaio = f f d^r, 

where {p)haio is the halo average density, and [Bhaio] = 
jBhaio^MdMlJ^MdM. 

As we can see, the boosting is determined by three 
functions: the halo mass function p{M, z), the boost for 
a single halo Bhaio as a function of the halo concentration 
Cvir, and the concentration dependence Cvir(M, z) of halo 
mass and z. 

TO _ _ ~ 1/(1 + Cyir)^) _ , , 

" 9(ln(I + Cvir) - Cvir/il + CyirW ^ ’ 

for the case of NFW halos. 
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FIG. 3: The annihilation signal J from the Milky Way as a 
fnnction of the angle 6 between the observation line and the 
direction towards the Galaxy center. The normalization of J 
is arbitrary. The solid line corresponds to the case when the 
clnmpy structnre disappears at rdestr = 50 kpc, the dashed 
line - to rdestr = 10 kpc. The boost factor at 50 kpc is set to 
be equal to Gso = 100. 


B. The boosting in the case of cored clumps 

What changes in the picture, if we consider the above- 
stated model instead of the traditional NFW one? It 
is surprising that the mass function should be almost 
the same for the cuspy and cored clumps. Indeed, two 
main factors defining p{M, z) are the spectrum of the ini¬ 
tial cosmological perturbations and the clump destruc¬ 
tion during the hierarchical structure formation. The 
first factor has nothing to do with the profiles of formed 
clumps. The second process (in the framework of the 
standard ACDM cosmology) is also not sensitive to the 
clump density profile. The objects of different masses col¬ 
lapse almost simultaneously in ACDM, and, as we have 
already mentioned, the hierarchical destruction time is 
shorter than the virilization time. Therefore, the ap¬ 
proach used by [5] and considering non-virialized clumps 
is exactly valid for in our case and leads to the same mass 
function. 

The only mechanism that significantly changes 
dp{M, z)/dM and essentially depends on the clump den¬ 
sity profile is the astrophysical clump destruction. Let us 
consider the Local Group. As we could see, the clumps 
situated farther than 50 kpc from the Milky Way center 
should survive under the most pessimistic assumptions 
about their structure. However, ^ 70% of the Milky 
Way dark matter lies outside of this radius [31], and we 
may suppose that the same situation takes place in the 
Andromeda galaxy. Moreover, along with the DM haloes 
of large and dwarf galaxies, the Local Group may contain 
an approximately equal quantity of dark matter that is 
not bound in the galaxies and presumably forms a large 
envelope [11, 21, 29]. Thus, ~ 85% of the Local Group ut¬ 
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FIG. 4: The annihilation signal J from the Milky Way as a 
function of the angle 9 between the observation line and the 
direction towards the Galaxy center. The normalization of J 
is arbitrary. The solid line corresponds to the case when the 
clumpy structure disappears at Vdestr = 50 kpc, the dashed 
line - to Vdestr = 10 kpc. The boost factor at 50 kpc is set to 
be equal to Gso = lO"*. 


terly kept its clumpy structure, and if we consider a large 
area of dark matter, significantly exceeding the maximum 
halo size, the modification of the averaged dp{M, z)/dM 
is not very significant. 

Of course, B^aio depends on the clump density profile. 
However, even in the case of the NFW clump with infinite 
central density, the annihilation signal (oc p^ ■ oc 

const) has no peculiarity in the center. Since the clear 
discrepancy between cores and cusps occurs only at r <C 
Ts, Bhaio for these cases may differ only on a factor < 
2 (see minimal model in [25] as an illustration). Thus 
the difference in value of Bhaio between cored and cuspy 
clumps is much less then the uncertainties of Cmr{M, z). 

Indeed, the uncertainty of Cvir{M,z) (especially the 
dependence on M) is the weakest point of any DM in¬ 
direct search theory. Since dp{M, z)/dM rapidly drops 
with M, the preponderant contribution into [Bhaio] is 
given by the smallest clumps, and we need to know their 
concentration. Typically Cmr{M,z) is defined from the 
results of N-body simulations that suggest a power-law 
dependence [30]. However, the ratio of masses of the 
smallest halos with properly resolved density profiles to 
the largest ones hardly exceeds 3 — 4 orders of magnitude 
even in recent large simulations containing ~ 10® test 
bodies. The mass hierarchy of the real Universe covers at 
least 21 orders of magnitude. Therefore, the estimations 
of Cair for the smallest clumps are based on a far interpo¬ 
lation. It is questionable whether the power-law behavior 
of Cair {M, z) holds within such a vast range of masses [4]. 
Recent N-body simulations suggest that Cmr < 70 even 
for the smallest clumps [13]. Moreover, Bhaio cx if 
Cvir ^ 1, i-e., Bhaio is sensitive even to small changes of 
Cair- All these facts make the predictions of [Bhaio] very 
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dependent on the suppositions of the Cmr{M,z) behav¬ 
ior. Now we will discuss several possible scenarios by the 
example of the Milky Way. 

Thus, if we consider a large area of dark matter, signif¬ 
icantly exceeding the maximum halo size, the difference 
between cuspy and cored scenarios can be surprisingly 
small, especially if we take into account the large un¬ 
avoidable uncertainties. For instance, if we could mea¬ 
sure the extragalactic background produced by the DM 
annihilation, we would be unable to choose between the 
scenarios because of our ignorance about Cyir{M, z) of the 
small clumps. However, if we consider a separate halo, 
the predictions of the models can be drastically different. 

In order to illustrate this, we consider the Milky Way 
halo, as it is a popular object for the indirect DM search. 
For simplicity, we assume that the clumpy structure ap¬ 
pears stepwise at some radius Tdestr- there is no substruc¬ 
ture inside this radius, and the substructures utterly sur¬ 
vive outside of Vdestr- In order to take into account the 
uncertainties of Cmr determination, we presume the total 
boost factor C 50 at r = 50 kpc from the Galaxy center. 
The absolute value of the signal depends on the particle 
nature of the DM, and we choose the normalization arbi¬ 
trarily (but the same for all the cases). The figures (2-4) 
show the annihilation signal J from the Milky Way as a 
function of the angle 6 between the observation line and 
the direction towards the Galaxy center. The solid line 
corresponds to Tdestr = 50 kpc, which is roughly consis¬ 
tent with the scenario of the least stable cored clumps, 
the dashed line — to Tdestr = 10 kpc, which approxi¬ 
mately models the NFW clumps [7]. Fig. 2, 3, 4 corre¬ 
spond to C 50 = 20 , 100 , 10 ^, respectively. 

As we can see, the difference between the scenarios is 


not significant, if the boosting is not very strong (C50 = 
20, Fig. 2). The dashed line on Fig. 3 approximately co¬ 
incides with Fig. 5 in [7]: the model with Tdestr = 10 kpc 
and C 50 = 100 indeed corresponds to that considered in 
this paper and assuming the NFW clumps and holding 
the power-law behavior of Cmr{M) obtained in N-body 
simulations. The cases of C 50 = 20 and C 50 = 10^ de¬ 
scribe slower and faster grows of Cmr{M), respectively. 
One can see that the difference between the cuspy and 
cored clumps becomes quite significant if Cso > 100 : the 
external regions of the halo generate a bright background, 
and it is significantly dimer if the clumps are easy to be 
destroyed. The background significantly modihes the an¬ 
gular profile of the signal. However, even if the boost fac¬ 
tor is as large as = 10 ^, the background produced by 
the outer clumpy regions of the halo is not bright enough 
to overshadow the Galaxy center. 
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